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摘要 ”本文 主要 介绍 了 近年 来 泛 函 重 整 化 群 (Functional Renormalization Group, fRG) 在 量子 色 动 力学 
(Quantum Chromodynamics, QCD) 相 变 及 其 临界 现象 研究 中 的 一 些 进 展 ,主要 侧重 于 与 QCD 的 临界 终点 
(Critical End Point, CEP) 相 关 的 非 微 扰 临 界 指数 的 计算 、 重 子 数 的 涨 落 等 。fRG 是 一 种 非 微 扰 连 续 场 论 的 理论 
方法 ,其 基本 思想 是 通过 重 整 化 群 能 标 从 高 能 向 低能 的 演化 ,连续 不 断 地 将 不 同 尺度 的 量子 涨 落 以 及 热 涨 落 积 
分 到 理论 体系 当中 。 本 文 讨论 了 非 微 扰 有 效 势 重 整 化 群 流 方程 和 不 动 点 方程 的 各 种 求解 方法 ,包括 泰勒 展开 、 
空间 维 数 e = 4 - 4d 展开 以 及 最 近 被 提出 的 非 局 域 势 直接 求解 方法 。 还 系统 地 讨论 了 与 QCD 的 临界 终点 等 临 
界 现 象 密切 相关 的 重子 数 的 涨 落 ,探讨 了 实验 观测 到 的 净 质 子 数 分 布 的 峰 度 ,也 就 是 四 阶 矩 随 碰 撞 能 量 非 单调 
依赖 关系 的 可 能 原因 。 
关键 词 QCD 相 结构 ，QCD 临界 终点 ， 泛 函 重 整 化 群 ， 临 界 现象 ， 相 对 论 重 离子 碰撞 
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Abstract Recent progress in studies on quantum chromodynamics (QCD) phase transition and related critical 
phenomena within the functional renormalization group (fRG) approach were reviewed, including the 
nonperturbative critical exponents and baryon number fluctuations, which are pertinent to the critical end point (CEP) 
in the QCD phase diagram. The fRG is a nonperturbative continuum field approach, in which quantum thermal 
fluctuations are successively integrated with the evolution of the renormalization group (RG) scale. Different 
methods of finding solutions to the flow or fixed-point equations of a nonperturbative effective potential have been 
discussed, for example, the Taylor expansion, expansion of the spatial dimension €= 4 — d, and the recently 
proposed direct solution of the global potential. Furthermore, the baryon number of fluctuations is relevant to the 
critical phenomena of the CEP. Both have been discussed, and one explores the underlying reasons for the observed 
non-monotonic dependence of the kurtosis of the net proton number of distributions on collision energy in 
experiments. 

Key words QCD phase structure, QCD critical end point, Functional renormalization group, Critical phenomena, 


Relativistic heavy ion collisions 


国家 自然 科学 基金 面 上 项 目 (No.12175030) 资 助 

第 一 作者 : 尹 诗 ， 男 ，1996 年 出 生 ，2018 年 毕业 于 大 连理 工大 学 ， 博 士 研究 生 ， 研 究 领 域 为 中 高 能 核 物理 
通信 作者 : 付 伟 杰 ，E-mail: wjfu@dlut.edu.cn 

收 稿 日 期 : 2022-12-22， 修 回 日 期 : 2023-02-22 

Supported by National Natural Science Foundation of China (No.12175030) 


五 


First author: YIN Shi, male, born in 1996, graduated from Dalian University of Technology in 2018, doctoral student, focusing on medium and high 
energy nuclear physics 

Corresponding author: FU Weijie, E-mail: wjfu@dlut.edu.cn 

Received date: 2022-12-22, revised date: 2023-02-22 


040002-1 


202306.00144v1 


chinaXiv 


核 技 术 


ChinaXiv 合 作 期 刊 


2023, 46: 040002 


近年 来 ,高 温 高 密 等 极端 条 件 下 强 相 互 作用 物 
质 的 性 质 、 量 子 色 动力 学 (Quantum 
Chromodynamics, QCD) 相 图 和 相 结 构 已 成 为 中 高 
能 核 物理 、 相对论 重 离子 碰撞 、 致 密 星 体 以 及 引力 波 
等 领域 的 前 沿 研究 热点 中 。 图 1 是 一 张 比较 典型 的 
QCD 相 图 的 示意 图 , 纵 坐 标 是 温度 , 横 坐 标 是 重子 
化 学 势 ,化 学 势 越 大 ,QCD 物质 的 密度 也 越 高 。 由 
图 1 可 以 看 到 ,在 高 温 低 密 的 区 域 相 边 界 是 一 条 虚 
线 , 这 条 虚线 不 是 一 条 严格 意义 上 的 相 变 线 , 而 是 一 
种 连续 的 变化 ,也 被 称 为 连续 过 渡 ,这 也 被 第 一 性 原 
理 的 格 点 QCD 模拟 所 证 实 中 ,同时 也 符合 实验 测量 
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图 1 QCD 相 图 的 示意 图 
Fig.1 Schematic phase diagram of QCD 

随 着 密度 或 者 化 学 势 的 增加 ,由 于 符号 问题 , 格 

点 QCD 的 计算 变 得 越 来 越 困 难 ,其 可 靠 性 通常 被 限 
制 在 图 1 中 js/7T< 2~3 的 区 域 " ,这 里 心 是 重子 化 
学 势 ,7 是 温度 。 其 他 理论 方法 ,比如 第 一 性 原理 的 
泛 函 连续 场 论 5、 低能 有 效 模 型 “的 计算 表明 ， 
在 高 密 区 域 连续 过 渡 有 可 能 转变 为 一 级 相 变 ,如 图 
1 中 红色 实 线 所 示 。 而 这 条 一 级 相 变 线 的 终点 ,也 
被 称 为 QCD 的 临界 终点 (Critical End Point, CEP) , 
也 就 是 图 1 中 黑 线 空心 小 圆圈 所 代表 的 点 ,将 一 级 
相 变 线 和 连续 过 渡 连 接 起 来 。 从 这 里 我 们 可 以 发 现 
CEP 在 QCD 相 图 上 的 特殊 地 位 , 它 的 存在 与 否决 定 
了 高 密 区 域 是 否 存 在 一 级 相 变 。 假 如 相 图 上 确实 存 
在 CEP ,那么 它 在 相 图 上 的 位 置 必然 能 够 反映 有 限 
温度 有 限 密度 非 微 扰 QCD 的 性 质 ,通过 实验 测量 这 
将 会 给 我 们 提供 一 个 研究 强 关 联 QCD 物质 性 质 的 
一 个 新 的 窗口 。 我 们 知道 ,作为 一 级 相 变 线 终点 的 
CEP 是 一 个 二 级 相 变 , 从 重 整 化 群 随 能 标 演化 的 角 
度 来 讲 ,CEP 所 对 应 的 系统 处 在 某 个 重 整 化 群 跑 动 
的 不 动 点 (也 称 为 固定 点 ) 的 临界 曲面 上 ”"。 我 们 在 
图 1 CEP 的 周边 示意 性 地 画 出 了 这 个 抽象 的 参数 空 
间 ,其 中 若干 黑色 带 箭 头 的 曲线 代表 参数 空间 中 重 
整 化 群 的 流动 ,可 以 看 到 CEP 确实 正好 处 在 重 整 化 


群 演化 的 不 动 点 上 。20 世 纪 70 年 代 初 ,自从 Wilson 
等 ”将 重 整 化 群 的 观点 和 方法 引入 到 二 级 相 变 和 
临界 现象 的 研究 中 ,相关 的 领域 得 到 了 突破 性 的 发 
展 。 现 在 我 们 知道 二 级 相 变 的 临界 行为 ,比如 刻画 
临界 行为 的 各 种 临界 指数 具有 普 适 性 ,其 只 依赖 于 
体系 的 空间 维 数 、 对 称 性 等 等 ,而 与 相互 作用 的 细节 
无 关 。 具 体 到 QCD 相 图 上 的 CEP, 它 作为 一 级 相 变 
线 的 终点 属于 Z(2) 对 称 性 普 适 类 ,也 就 是 伊 辛 模型 
的 普 适 类 。 
目前 ,有 很 多 重 离子 实验 设备 正在 或 者 准备 在 
实验 中 探寻 CEP™””””。 我 们 知道 ,在 三 级 相 变 点 附 
近 , 体 系 的 关联 长 度 明 显 增加 ,在 理想 热力 学 极限 的 
情况 下 ,CEP 处 的 关联 长 度 趋 于 发 散 。 热 力学 和 统 
计 物 理 的 基础 知识 告诉 我 们 ,有 一 些 观测 量 ( 如 各 种 
守恒 荷 的 涨 落 ? 对 这 种 临界 动力 学 所 导致 的 奇异 性 
特别 敏感 ,而 且 阶 数 越 高 涨 落 越 敏感 。 早 在 二 十 多 
年 前 ,就 提出 在 重 离子 碰撞 实验 中 可 以 利用 守恒 荷 
涨 落 随 碰撞 能 量 的 非 单调 依赖 关系 来 探寻 QCD 相 
图 上 的 CEP”” ,也 可 以 参考 文献 [26]。 过 去 10 年， 
鲁 克 海 文 国家 实验 室 相 对 论 重 离子 碰撞 实验 
(Relativistic Heavy Ion Collider, RHIC) 束 流 能 量 扫 
描 项 目 (Beam Energy Scan, BES) 第 一 阶段 的 实验 ， 
测量 了 净 质 子 、 净 电荷 以 及 净 KK 介 子 多 重度 分 布 的 
各 阶 矩 及 其 它们 之 间 的 关联 ” ,取得 了 丰硕 的 成 
果 , 特 别 是 在 最 近 的 实验 中 ,以 3.1lc 的 显赫 性 观测 到 
了 净 质 子 多 重度 分 布 的 峰 度 随 对 撞 能 量 的 非 单调 依 
HRR", 
理论 研究 方面 ,前 面 我 们 已 经 提 到 由 于 有 限 化 
学 势 情况 下 的 符号 问题 , 格 点 QCD 模拟 被 限制 在 
Ls/T7< 2~3 的 区 域 , 但 在 这 个 区 域 并 没有 发 现 
CEP'“”!。 近 年 来 ,第 一 性 原理 泛 函 连续 场 论 的 方法 
得 到 了 快速 的 发 展 " 中 , 泛 函 连续 场 论 的 计算 不 仅 
在 化 学 势 比较 小 的 区 域 ,例如 js/T < 2, 得 到 了 格 点 
QCD 的 检验 ,并 且 其 计算 的 可 靠 性 还 被 拓展 到 
Ha < 3~4 的 区 域 。 此 外 ,近年 来 不 同 的 连续 场 论 
WN Awe , mW yw eh Me RM 化 群 (Functional 
Renormalization Group, fRG) 以 及 Dyson-Schwinger 
方程 都 对 CEP 在 相 图 上 的 位 置 进行 了 估计 ,这 两 种 
方法 给 出 了 在 误差 范围 内 一 致 的 结果 ,都 表明 CEP 
所 对 应 的 化 学 势 约 在 600 MeV, 
本 文 计 划 介 绍 近年 来 泛 函 重 整 化 群 在 QCD 相 
变 及 其 临界 现象 研究 中 的 一 些 进展 ,包括 非 微 扰 临 
界 指数 的 计算 、 重 子 数 的 涨 落 等 。 保 G 是 一 种 非 微 
扰 连续 场 论 的 理论 方法 ,其 基本 的 思想 是 通过 重 整 
化 群 能 标 从 高 能 向 低能 地 跑 动 ,连续 不 断 地 将 不 同 
尺度 的 量子 涨 落 以 及 热 涨 落 积 分 到 我 们 的 体系 当 
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HAS) ARG 已 经 被 应 用 到 非 微 扰 的 很 多 领域 ,QCD 
相关 的 保 G 综述 可 参考 文献 [1$,$1-$8]。 近 年 来 的 
一 些 重要 进展 可 参考 文献 [12 ,3$9-65]。 


1 函 重 整 化 群 简介 


首先 ,生成 泛 函 如 下 : 
ZJ]= [(Do)exp| -STG6]- ASLE] + J°6,\1) 


式 中 :多 代表 了 某 一 场 论 中 所 有 的 场 ;S[$] 是 其 经 
典 作 用 量 。 在 这 里 引入 一 个 能 标 k 依 赖 的 红外 抑制 
函数 AS,[ 人 $], 其 作用 是 抑制 生成 泛 函 (1) 中 动量 模 
式 g < 的 量子 涨 落 ,而 不 影响 g > 的 涨 落 。 通 常 ， 
为 了 简便 选用 双 线 性 场 形式 的 红外 抑制 函数 : 
AS,[O]= 5 6,R6, (2) 
其 中 :对 于 玻 色 型 的 指标 有 R2% = R”, 对 于 费 米 型 
R? =- R”, 
下 面 以 只 有 一 个 分 量 的 标量 场 p 为 例 , 其 红外 
抑制 函数 为 : 


ASi[9]= 5 atxaty p(X) R(xy)o(y) (3) 


利用 动量 空间 的 传 里 叶 变换 
oy = foo (4) 

以 及 ; 
Rov- ER)" 加 

得 到 : 
Asto1= 3 So- Rg)o(g) © 
根据 前 面 的 要 求 , 希 望 抑制 函数 具有 如 下 的 渐 

近 行为 ， 


Risoli vo Rea (fixed g) (7) 

此 外 ,为 了 有 选择 性 地 抑制 低 动 量 区 域 的 量子 
涨 落 ,进一步 要 求 : 

Ri(g)| K, Ri(g) (8) 

显然 ,满足 上 面 要 求 的 抑制 函数 有 无 穷 多 种 选 
择 ,下 面 介绍 两 类 在 文献 中 使 用 比较 多 的 抑制 函数 。 
第 一 类 是 指数 型 的 抑制 函数 : 

REG) = Fro E) rho( = — 0) 

参数 "控制 抑制 函数 在 9 ~ 处 的 陡峭 程度 。 
另 一 类 更 适合 解析 计算 的 抑制 函数 被 称 为 平坦 抑 
函数 ,或 者 称 为 最 优化 抑制 函数 re 


ek 


= 


REQ) = Pral), ras) =[4- tod =x) 
(10) 
式 中 :0(x) 是 阶梯 函数 。 
由 场 论 的 一 般 理论 ,知道 连通 关联 函数 的 生成 
ZETANS R: 
W[J1]= InZl 7] (11) 


这 也 被 称 为 Schwinger 函数 , 场 的 期 待 值 可 以 通 


过 对 Schwinger 函数 求 外 源 的 一 阶 导 数 得 到 : 
= _ Wis) 
®,=(@,) a (12) 
再 求 一 次 导数 可 以 得 到 传播 子 : 
/SS a(S S Ow 
Gn = (5), = (6.4) -(4,)(0,) = PHL 
(13) 
对 Schwinger 函数 作 Legendre 变换 ,得 到 单 粒子 


不 可 约 (1P1) 


TAHE: 
TP] = WJ] tJ B,- AS [TP] 04 

引入 如 下 记号 : 

J B, = yD, J" (15) 


H 


( 


RLN FE: 


y 


内 中 : 当 a,b 是 费 米 型 指标 ,y; = 0%; a, DERE 
型 指标 ,ys = 5654。 这 样 得 到 : 


6(T[B]+ AS[G]) 


= yJ’ (16) 
进而 可 以 得 到 传播 子 和 有 效 作用 量 的 关系 : 
Ge=X(TPIG]+ASPIG]) (7) 


这 里 有 : 


FP[®]+ ASP[#))" = 


o(I,[]+ AS,[®]) 
ôb, ôD, ue) 
由 式 (1) 和 (11) 容 易 得 到 Schwinger 函数 随 大 的 


aW,LI1== 5Sr[(0,R,)G.] ~ 32RD, (19) 
这 里 引入 了 重 整 化 群 时 间 1= In (KA),4 是 某 


一 参考 能 标 , 比 如 紫外 演化 初始 能 标 。 在 式 (19) 中 
还 引入 了 一 个 记号 超 求 迹 : 
STr[(a,R,)G; | = (RE) Grea (20) 
再 一 次 利用 Legendre 变换 : 


0 


LLP] = 0W,LJ]- 


a,AS,[®]= Z STr[(3,R.)G] 


(21) 
也 被 称 为 


这 就 是 有 效 作 用 量 的 流 方程 ， 


Wetterich 方程 5 。 
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有 效 作用 量 的 流 方程 (21) 可 以 改写 为 : 
a0, ]= 5STr[i In(TOIG]+R)] (22) 
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其 中 :6, 代 表 只 对 抑制 函数 的 重 整 化 群 能 标的 依赖 
性 求 偏 导 ,这 a aa 


(TPs)? = Fo Als (23) 


这 里 使 用 了 左右 导数 ， a 
所 多 出 来 的 负 号 。 作 如 下 的 组 合 : 
1 


TO[G]+R= P+ F (24) 
其 中 :了 是 两 点 关联 函数 矩阵 ,不 含有 相互 作用 ,其 
到 数 就 是 传播 子 矩 阵 ;三 是 相互 作用 和 矩阵 ,也 就 是 相 
互 作用 项 点。 将 上 面 有 效 作用 量 的 流 方程 作 形式 上 
的 泰勒 展开 : 


ed 


1 


a0,= 5ST r[ð, nee 7 STrd,InP + > ;small Fr 本 
4 (25) 
psi oF + stra 1 F) - str 5 F +e 
it He E E 7 TE HR AID EA, 1 方程 的 右边 需要 计算 动量 空间 中 的 两 点 函数 ; 
到 各 阶 关联 画 数 ,包括 各 种 传播 子 和 顶点 的 流 方程 。 row(gigj= FTI] [zg + 7(p))8, + 
事实 上 ,还 有 一 种 更 简单 的 方法 可 以 很 方便 地 plaga) ~ Fee" 
给 出 各 阶 关 联 函 数 的 流 方程 ,可 以 借用 同 微 扰 论 完 V(p)6.0,| (2n)'5"(q' + q) (30) 


全 一 致 的 费 曼 规 则 。 首 先 我 们 定义 一 个 一 般 的 n 点 
1PI 关 联 函数 ,或 者 称 为 n 点 顶点 : 

oT LP] 
df, 5D, 


Veo ti, 三 一 an = | (20) 


P=< p> 
这 样 顶 点 ”的 流 方程 可 以 表示 为 : 


n | all one - loop correction | 
7) 


OV, 


kb, 0” 


=0 . 
| diagrams of V7)... 


更 详细 的 讨论 可 参考 文献 [15]。 
2 O(N ) 标 量 理论 


在 连续 相 变 的 临界 点 附近 ,系统 的 关联 长 度 会 
趋 于 发 散 ,系统 的 微观 细节 会 被 临界 点 附近 强烈 的 
涨 落 掩盖 ,临界 点 附近 物理 量 的 行为 通常 可 以 用 临 
界 指数 来 描述 ,它们 只 依赖 于 系统 的 空间 维 数 与 对 


选择 场 的 零 分 量 ( 纵 向 ) 来 携带 场 的 期 待 值 而 让 
其 他 方向 (横向 ) 的 期 待 值 为 零 。 
mix = 20V(p)dn + Vilp) 
Mon = Mi-o = 2VE(p) + Vip) (31) 
Mis = Mso = Vi(p) 
动量 空间 中 的 抑制 函数 写 为 : 
Rg (a'g) = Ri(q) Ore (q + q)o, 
=Zyqr(q'/k?)(2n)'6(q'+q)d; (32) 
抑制 函数 对 重 整 化 能 标的 导数 为 : 
,RI(q) = Zorxk (nmr (x) + 2xr'"(x)) (33) 
其 中 ,x = gik. REEN: 
_ Olyr 
Lok 
把 这 些 代 入 Wetterich 方程 并 做 如 下 无 量 纲 化 


(34) 


称 性 。d 维 的 O(N ) 理 论 则 恰好 涵盖 了 不 同 的 空 
维 数 与 对 称 性 ,在 临界 点 附近 它 可 以 描述 大 部 分 
实物 理 系统 的 二 阶 相 变 。 
Wilson 的 重 整 化 群 * 则 是 描述 相 变 的 一 个 强 
大 的 理论 工具 。 二 阶 相 变 由 重 整 化 群 流 方程 所 确定 
的 临界 面 上 的 稳定 固定 点 所 描述 。 ERG 则 更 进一步 
地 在 重 整 化 群 的 流 中 包括 了 非 微 扰 的 信息 ,这 让 我 
们 能 在 空间 维 数 比 较 小 时 更 精确 地 描述 相 变 。 
首先 来 看 O(N ) 标 量 理 论 的 有 效 作用 量 : 


Plo] = fax E Za (V$) + Vip) - nh (28) 


式 中 :h 是 外 场 ;Vi(p) 是 一 个 关于 场 的 一 般 函 数 , 通 
常 也 被 称 为 有 效 势 ,其 中 有 


p= #2, p = ay (29) 
显然 ,Vi(p) 具 有 O(N) 对 称 性 。 


在 Wetterich 


重 定义 : 
p= EW Oe u,(p) = EVP), 
-2 _ mM m= max (35) 
2 ”Zk 
可 以 得 到 无 量 纲 化 的 有 效 势 的 流 方 程 : 
du(p) = du(p)+(d-2+ n)pu'(p) + 
C [1(m2) + (N - 1) (m2) ] (36) 
这 里 函数 1 的 具体 形式 与 抑制 函数 的 选择 有 
关 。 原 则 上 来 说 物理 量 在 固定 点 附近 的 行为 不 依赖 
于 抑制 函数 的 选取 ,但 是 ,在 实际 的 计算 中 通常 要 对 
我 们 的 理论 做 截断 , 璧 如 通常 用 的 局 域 势 近似 
(Local Potential (ee LPA) 近 似 就 是 只 考 
虑 梯度 展开 的 最 低 阶 ,这 会 引入 一 定 程度 上 的 抑制 
函数 的 依赖 性 。 这 时 要 选择 一 个 恰当 的 抑制 函 
a 论 尽 量 减少 这 种 依赖 性 ,在 低 阶 截断 的 情况 
也 能 很 好 地 描述 真实 的 物理 。 关 于 抑制 函数 的 选 


040002-4 


202306.00144v1 


chinaXiv 


P 诗 等 ; 临界 现象 与 泛 函 重 整 化 群 


ChinaXiv 合 作 期 刊 


择 及 其 优化 ,更 详细 的 讨论 可 参考 文献 L66-67]。 在 
这 里 采用 式 (10) 中 介 ee 
r(x)= (x =e -x) (37) 
此 时 有 I(x)= 1(1 ea 
1 1 1 2 
~ 2 (4ny® TOA) (2 anla Ty Zle» 


3 ”高 斯 固定 点 


暂时 不 考虑 反常 量 纲 , 设 7 = 0, 并 将 方程 (36) 
在 万 = 0 附近 作 泰 勒 展 开 ” : 


== An TR 
UP) È F (39) 


其 中 :展开 系数 ,对 应 着 没有 动量 依赖 的 无 量 纲 化 
的 2n 点 耦合 强度 ,特别 地 ,1 是 场 为 零 时 的 曲率 质 
量 。 如 果 展 开 到 二 阶 并 且 忽 略 四 点 相互 作用 对 四 点 
强度 1 的 流 的 贡献 ,得 到 
aA, = 24,-C(N+2)A,(1+4,) 
ðA, = (4 - d) h 
这 样 就 给 出 了 相应 耦合 强度 的 8 函数 : 
0,4, = B,(A) (41) 
SB = B,= 0, ETRE Fh aD 77 Fe ate T 
个 代数 方程 , 它 的 解 是 重 整 化 群 方程 的 固定 点 : 
=A =0 (42) 
如 果 重 整 化 群 方程 的 固定 点 满足 4; = 0, 我 们 
说 这 样 的 固定 点 是 平庸 的 或 者 说 是 高 斯 固定 点 。 一 
般 来 说 , 像 方程 (40) 这 样 的 重 整 化 群 流 方程 不 能 解 
析 求 解 ,在 固定 点 附近 利用 小 量 : 
= 一， R= {hda} (43) 
可 以 将 方程 (40) 线 性 化 : 


4 > > (=o CWE 
0,67 = M67 = 
0 -(4-d) 


(40) 


(a) 9 


(a) 


Critical Line 


—0.02 —0.01 0 0.01 0.02 
4, 


图 2 高 斯 固定 点 附近 重 整 化 群 


XEEN T ERE: 
m, = 2È (45) 


MERR i 4S AS IE (EPP AS iE Bt} AI y, 
和 e, 上 面 重 整 化 群 流 方程 的 解 可 以 一 般 地 写成 : 
Oh = Dae, (46) 


式 中 :w; 是 常数 小 量 。 对 于 这 里 讨论 的 高 斯 固定 点 ， 
它 的 本 征 值 和 本 征 矢 量 分 别 为 ; 
Yı = 2,y,=- (4 - d) 
é,=(1,0)".2, -(- 1 
这 样 可 以 得 到 与 关联 长 度 相 关 的 临界 指数 v = 
-1/y1= 0.5。 图 2 展示 了 重 整 化 群 方程 (40) 的 流 在 
和 =-4 罗 平面 上 的 演化 行为 ,箭头 所 指 的 方向 是 重 整 
化 群 向 红外 流动 的 方向 ,yi 是 唯一 的 一 个 高 斯 回 定 
点 ,el, 2; 是 两 个 本 征 矢量 。 由 于 重 整 化 的 能 标 向 红 
外 流动 , 即 式 (46) 中 上 < 0, 所 以 , 当 本 征 值 大 于 零 
所 对 应 的 本 征 矢 量 流向 固定 点 ,而 当 本 征 值 小 于 零 ， 
所 对 应 的 本 征 矢 量 流出 固定 点 。 图 2(a) 是 空间 维 数 
d > 4 时 的 流 图 ,在 这 种 情况 下 y, < 0,y, > 0, 因 而 可 
以 看 到 , 重 整 化 群 的 流 沿 本 征 矢量 2 流出 高 斯 固定 
点 ,而 沿 2, 流 入 固定 点 ,这 样式 (46) 中 相应 的 参数 a， 
是 相关 参数 ,a, 是 无 关 参 数 。 我 们 称 像 这 样 只 有 一 
个 相关 参数 的 固定 点 为 稳定 固定 点 , 它 对 应 于 物理 
上 的 三 阶 相 变 。 参 数 空间 中 相关 参数 为 零 的 子 空间 
也 被 称 为 临界 面 , 或 者 临界 线 , 显 然 图 2(a) 中 本 征 矢 
量 人 2, 即 为 临界 线 。 随 着 重 整 化 群 流 的 跑 动 ,无论 从 
哪个 点 出 发 ,无 关 参 数 的 值 会 趋 于 零 ,唯一 的 相关 参 
数 的 值 会 发 散 。 由 此 得 出 结论 ,当空 间 维 数 d > 4 
时 ,高 斯 固定 点 是 唯一 的 固定 点 ,也 是 物理 的 固定 
点 ,dg=4 是 上 临界 维 数 , 上 临界 维 数 之 上 ,临界 指数 


" (47) 


0.0 - < 
-0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 


方程 的 流 图 (a)d>4,(b)2<d<4 


Fig.2 Flow diagram of RG equations around the Gaussian fixed point in dimensions (a) d > 4, (b)2 < d < 4 
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与 平均 场 理论 相同 。 图 2(b) 是 空间 维 数 2< dg< 4 
时 的 流 图 ,可 以 看 到 ,本 征 矢量 2, 和 e, 都 是 从 固定 点 
流出 ,此 时 的 高 斯 固定 点 有 两 个 相关 参量 ,是 不 稳定 
的 固定 点 ,因而 2 < d < 4 情况 下 的 高 斯 固定 点 不 是 
物理 的 固定 点 ,需要 寻找 新 的 固定 点 。 

E d < 4 时 无 量 纲 化 之 前 的 四 点 相互 作用 强度 
也 是 相关 参数 , 它 不 能 被 忽略 。 更 一 般 地 ,27 顶点 
耦合 强度 (无 量 纲 化 之 前 ) 的 量 纲 为 [入 ]= 2n - (n - 
1)4, 通 过 假设 ,是 边缘 参数 , 即 [7 ,1]= 0, 可 以 
aa e ee d,=2+2/n,4d,<d<d,,, 
时 有 n 个 相关 参数 ,更 完整 的 计算 表明 ,此 时 在 高 斯 
EE EUn AT EREHE, Aen RTN 固定 
点 的 本 征 矢量 。 特 别 地 , 当 g= 2 时 ,在 高 斯 固定 点 
旁边 的 本 征 扰动 有 无 穷 多 个 小 于 零 的 本 征 值 ,因而 
有 无 穷 多 个 相关 参数 ,以 及 无 穷 多 条 流出 固定 点 的 
AGERE o 


4 8=4-4d 展 开 及 Wilson-Fisher 固 定点 


从 83 分 析 知 道 , 在 2 <&< 4 时 ,高 斯 固定 点 不 
能 描述 物理 上 的 二 阶 相 变 。 而 此 时 四 点 相互 作用 强 
已 经 是 相关 参数 ,因此 ,需要 考虑 它 对 式 (40) 中 
0,4, 的 贡献 ,这 样 方 程 (40) 变 为 : 
ðA, == 24, - C(N+2)A,(1+4,) 


0h = (4 - d)A, + 2C (N + 8) (A;) (1 + 4,) 


_ (48) 


这 时 M 和 矩阵 为 : 
M= 
N+2 N+2 
-2+ see | 
0 E 


(51) 
这 个 矩阵 的 一 个 本 征 值 及 其 对 应 的 本 征 矢 量 分 
别 为 : 


N+2 2 r 
y= 24+ 5s, é=(1,0) (52) 
yi < 0, 因此 这 是 相关 参数 ,对 应 的 临界 指数 y 
精确 到 Cle) 阶 为 "3(- 1/y 泰勒 展开 到 e 阶 ): 
vt (53) 
另 一 个 本 征 值 及 其 对 应 的 本 征 矢 量 为 y,= 
ee, =(- B, 1), Ry, > 0 是 无 关 参 数 , 这 里 的 常 


数 B 是 : 
N+2 
TIN zeh ¥ 
3 


1 N-4 
aa ~ bewen +n 
(54) 
图 3(a) 给 出 了 2 <d <4 情 况 下 重 整 化 群 方程 
(49) 的 流 图 ,可 以 清楚 地 看 到 ,在 平面 上 有 两 个 固定 
点 ,一 个 是 处 于 原点 位 置 的 高 斯 固定 点 ,显然 高 斯 固 


B= 5C(N+2)]I 


ETRY it HERR TA re OB XE ah A; = 2, = 0 之 
外 ,从 4 的 流 方程 ,可 以 发 现 , 当 e=4-qd<1, 有 一 
Aime, re 的 固定 点 ,在 这 个 固定 点 上 可 以 做 近 
似 1 +A; ~ 1 其 满足 的 方程 为 : 
ðA, =- 24; -C (N+ 2), =0 


49 
0,4, = el} + 2C (N + 8) (23) =0 (19) 
该 代数 方程 的 解 为 : 
»~_ N+2 pe E 
A= Gnas?” aea 


可 以 看 到 ,它们 都 是 Cs) 阶 ,这 个 非 零 的 固定 
点 被 称 为 Wilson-FisherC(W-F) 固 定点 。 在 W-F 固定 
点 附近 将 流 方程 (49) 线 性 化 ,并 且 严 格 计算 到 Oe) 
阶 和 


定点 是 不 稳定 的 固定 点 ,这 同 我 们 前 面 的 分 析 是 一 
致 ;另外 一 个 标记 为 巡 的 固定 点 是 WE 固定 点 , 它 是 
一 个 稳定 的 固定 点 ,对 应 物理 的 二 级 相 变 。 图 3Cb) 
在 对 称 群 参数 N 和 空间 维 数 4 构 成 的 平面 上 给 出 了 
W-F 固 定点 精确 到 O(e) 阶 的 临界 指数 y 的 计算 数 
值 , 即 式 (53)。 可 以 看 到 v 随 着 a 的 减 小 和 NV 的 增 大 
而 单调 增加 。 在 大 NN 极限 时, 球 模型 给 出 的 严格 结 
果 为 v= 1/(d -2), 由 于 方程 (53) 的 计算 结果 只 精 
确 到 O(e) 阶 ,在 d 比 较 小 时 , 它 显然 不 能 与 大 N 极 
限时 的 结果 吻合 。 

EON REEM y = 0, 可 以 利用 标 度 率 ， 
计算 出 准确 到 O(c) Gr HI Elle FES 


amami TESS a 
N+2 

=(2- = ] 十 

I 


5 mA AAT 


我 们 计算 了 精确 到 C(e) 阶 的 临界 指数 ,但 是 它 
会 在 维 数 比 较 小 时 失效 ,一 个 最 直接 的 改进 是 不 对 


E 3 
3v(d 21+) 2 2(N+8)° 本 
BPH 25 is 
d-2+ 


v= 1/y1 做 ge 展开 ,而 直接 保留 方程 (52) 中 yy 的 表达 
式 , 那 么 可 以 得 到 ; 
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v 
0.967 54 


图 3 (a)2< d< 4 情况 下 W-F 固 定点 附近 重 整 化 群 方程 的 流 


区 


b) 在 对 称 群 参 数 N 和 空间 维 数 4 的 平面 上 得 出 的 W-F 固定 点 


,(b 
精确 到 O(s) 阶 的 临界 指数 * 的 计算 数值 


Fig.3 (a) Flow diagram of RG equations around the W-F fixed point in 2 < d < 4, (b) critical exponent v in the order of O (e ) at the 
W-F fixed point obtained in the plane of N of the O(N ) symmetry and the dimension d 


E 8+N 
8-2N+d(2+N) 
在 取 N 一 % 极 限时 ,得 到 球 模型 的 严格 结果 
v=1/(d-2)。 根 据 MWH (Mermin-Wagner- 
Hohenberg) 定 理 '”"", 在 二 维 时 ,不 存在 连续 的 自发 
对 称 性 破 缺 ,因此 在 O(N ) 模 型 中 , 当 4d = 2,N>2 
时 ,系统 不 存在 二 级 相 变 ,临界 指数 vy 会 发 散 , 而 这 
里 计算 出 的 v 显 然 只 有 在 N 一 % 才 会 发 散 , 因 此 ,这 
里 计算 的 v 在 4 和 NV 比较 小 时 依然 存在 比较 大 的 
误差 。 
由 朴素 的 量 纲 分 析 可 知 , 无 关 参 数 要 求 
[4,] < 0, 结 合 之 前 对 2z 顶点 耦合 强度 的 量 纲 分 析 ， 
1X Db BER n> d/(d-2), 4 d=4,n>2;d= 
3,n>3;d=25,n>5;d=2.1,n>21;3°:-. © 
然 , 当 空间 维 数 越 低 ,就 需要 在 重 整 化 群 的 流 方程 中 
包含 越 来 越 多 的 耦合 常数 ,或 者 说 式 (39) 中 的 有 效 
势 的 泰勒 展开 阶 数 需要 越 来 越 高 。 这 里 计算 W-F 固 
定点 处 临界 指数 的 方法 是 利用 式 (39) 将 无 量 纲 化 的 
有 效 势 在 5 = 0 处 泰勒 展开 到 nn. 阶 ,然后 把 式 (39) 
代入 有 效 势 的 流 方程 (36) ,这 样 就 得 到 了 mm .个 浓 
微分 方程 ,进一步 令 9,4, = 0, 最 后 得 到 ni 个 代数 


(56) 


方程 。 从 这 些 耦 合 的 代数 方程 中 可 以 得 到 一 系列 固 
定点 ,并 且 在 每 一 个 固定 点 上 线性 化 这 些 常 微分 方 


FE, A APS S77 FE (44) PY MERE Th SE AE AY AS 
征 值 ,W-F [Al E OT VAY MERE HAA AS ASME 
值 小 于 零 ,这 一 个 小 于 零 的 本 征 值 的 绝对 值 就 是 我 


图 4 对 称 群 参数 一 % 极 限 下 , 取 N = 5000, 不 同 展开 阶 
数 计算 出 的 临界 指数 y 随 维 数 4 的 变化 
Fig.4 Critical exponent v in the limit of large N 一 oo as a 
function of the spatial dimension d calculated with different 
values of the expansition order ne 


在 图 5 中 ,计算 了 不 同 展开 阶 数 下 临界 指数 v 在 
d- 平面 上 的 行为 。 随 展开 阶 数 的 增加 , 临 界 指数 
确实 可 以 被 计算 到 越 来 越 小 的 维 数 ,但 是 不 能 任意 
地 接近 d = 2, 当 展开 阶 数 越 来 越 高 ,在 维 数 比较 小 
时 ,临界 指数 的 行为 越 来 越 奇异 。 这 是 由 于 有 效 势 
泰勒 展开 存在 收敛 半径 。 随 着 维 数 越 来 越 接近 2， 
有 效 势 最 小 值 点 的 位 置 会 越 来 越 远 离 零 点 ,而 由 于 
泰勒 展开 收敛 半径 的 问题 ,在 某 个 维 数 以 下 ,即使 我 
们 展开 到 无 穷 阶 也 不 可 能 把 有 效 势 展开 到 最 小 值 点 
以 外 。 想 要 解决 这 个 问题 ,我 们 可 以 利用 最 近 提 出 
的 一 种 方法 ,直接 求解 非 局 域 的 有 效 势 的 固定 点 


们 所 要 计算 的 临界 指数 v 的 倒数 。 图 4 展示 了 在 大 
N 极 限 下 ,不 同 展开 阶 数 计算 出 的 临界 指数 v 的 数值 
结果 及 其 对 空间 维 数 的 依赖 性 ,可 以 看 到 随 展 开 阶 
数 的 提高 ,临界 指数 的 结果 越 来 越 接近 真实 值 。 但 
是 不 论 展开 到 多 少 阶 临界 指数 随 4 的 变化 都 存在 非 
单调 性 行为 。 


方程 ,然后 在 固定 点 有 效 势 解 的 基础 上 利用 本 征 微 
扰 的 方法 , 求 得 相应 的 临界 指数 ,相关 的 更 多 细节 可 
参考 文献 [72]。 图 6 利用 这 一 方法 计算 得 到 的 临界 
指数 v, 相 较 于 图 4 的 计算 结果 可 以 看 到 , 非 局 域 求 
解 得 到 的 临界 指数 在 维 数 比 较 小 的 时 候 不 再 出 现 奇 
异 的 行为 。 利 用 这 一 方法 求解 得 到 的 临界 指数 只 存 
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在 由 于 截断 带 来 的 系统 误差 。 由 于 这 里 的 计算 没有 
考虑 反常 量 纲 的 影响 ,在 4= 2 时 ,临界 指数 vy 在任 
意 和 数值 下 都 会 发 散 , 这 比 MWH 定理 的 限制 更 强 。 
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如 果 进 一 步 考 虑 反常 量 纲 的 影响 ,临界 指数 数值 计 
算 的 系统 误差 会 进一步 降低 ,并 且 与 MWH 定理 的 
要 求 完全 符合 ,相关 的 研究 正在 进行 中 。 
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S 不 同 展开 阶 数 情况 下 计算 得 到 的 临界 指数 * 随 空间 维 数 4 和 对 称 群 参数 的 变化 


Fig.5 2D plot of the critical exponent v as functions of the spatial dimension d and the symmetry N calculated with different values 


of the expansition order n 


LPA-Exact 


4.0552 
3.320 1 


oo 
T 


2.7183 

2.225 5 
= 1.8221 
— 1.4918 
— 1.2214 
= 1.0000 
0.8187 


an 
T 


上 
T 


N 
T 


0.670 3 


6 LPA 近似 下 ,临界 指数 vy 的 精确 结果 及 其 随 空间 维 数 d 
和 对 称 群 参数 X 的 依赖 行为 
Fig.6 Exact results of the critical exponent v as functions of 
the spatial dimension d and the symmetry N with 
the truncation LPA 


6 BT RUNS 


我 们 讨论 在 相对 论 重 离 子 碰撞 实验 中 与 临界 现 
象 密切 相关 的 涨 落 测量 ,这 也 是 实验 中 寻找 QCD 临 


trunc 


界 终点 (CEP) 的 主要 手段 之 一 。QCD 的 临界 终点 是 
属于 Z(2) 对 称 性 普 适 类 的 二 级 相 变 ,也 就 是 伊 辛 普 
适 类 。 涨 落 随 关 联 长 度 的 究 次 随 涨 落 的 阶 次 增加 而 
增 大 ,因此 越 高 阶 的 涨 落 对 临界 现象 越 敏感 。 美 国 
布鲁克 海 文 国 家 实验 室 相 对 论 重 离子 碰撞 实验 
(RHIC) 束 流 能 量 扫描 项 目 (BES) 期 待 着 利用 守恒 
荷 ( 包 括 重 子 数 或 者 净 质 子 数 . 电 荷 以 及 奇异 数 等 ) 
的 高 阶 涨 落 及 其 对 碰撞 能 量 的 非 单 调 依赖 性 来 寻找 
QCD 相 图 上 的 CEP, 在 BES 第 一 阶段 的 实验 中 取得 
SRE BEN RR. 

为 解释 这 些 实验 数据 ,并 且 将 实验 测量 与 QCD 
相 变 以 及 CEP 联系 起 来 , 涨 落 的 理论 计算 与 预言 也 
必 不 可 少 。 近 年 来 , 涨 落 的 理论 研究 和 计算 也 取得 
了 很 多 重要 的 进展 ,譬如 关于 涨 落 的 格 点 QCD 计算 
可 参考 文献 L[73-79] , 保 G 框架 下 低能 有 效 模型 
(Low Energy Effective Theory,LEFT) 的 研究 可 参考 
文献 [80-92] ,相对 论 平均 场 可 参考 文献 [93-97] ， 
Dyson-Schwinger 方程 ”2%。 特 别 是 近年 来 在 保 G 
的 框架 将 低能 有 效 模型 与 第 一 性 原理 QCD 相 结合 ， 
构建 了 QCD 辅助 的 低能 有 效 理论 (QCD-assisted 
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LEFTs) ,并且 在 此 理论 框架 下 研究 了 重子 数 分 布 的 
(i 8 BEAU PE), 重子 -奇异 数 关联 , 超 高 阶 重 
子 数 的 涨 落 中 。 
6.1 QCD 辅助 的 低能 有 效 模型 下 的 重子 数 涨 落 
利用 QCD 辅助 的 低能 有 效 理论 计算 重子 数 的 
涨 落 ,并 且 将 计算 的 结果 与 实验 数据 进行 比较 。 
Polyakov 夸克 介子 (Polyakov-Quark-Meson, PQM) 
模型 非常 适合 用 来 计算 QCD 的 热力 学 量 。PQM 模 
型 中 包含 了 夸克 自由 度 以 及 6 介子 .x 介子 的 自由 
度 。 这 里 只 加 入 了 标量 和 硒 标 量 介子 ,因为 这 两 个 
介子 的 真空 质量 最 轻 , 在 描述 手 征 对 称 性 的 自发 破 
缺 与 恢复 中 起 到 了 决定 性 的 作用 ,其 他 质量 更 重 的 
介子 的 作用 都 可 以 视 为 次 领头 阶 的 贡献 。 在 夸克 介 
子 模型 的 基础 上 我 们 加 入 了 了 Polyakov 圈 的 贡献 ,以 
引入 夸克 禁闭 与 退 禁 闭 的 信息 。 这 样 就 可 以 给 出 低 
能 有 效 模型 的 有 效 作 用 量 : 
T= | {ZUP ~ yo + ig4o)]g + 


F 24(,8) + hg (tha + t: n)a + Vy Pedy) - 
co } (57) 
等 号 右边 第 一 部 分 为 夸克 的 运动 学 项 ,4 和 5 分 
别 代表 正 反 夸克 场 ,Z,, 为 专 克 的 波 函数 重 整 化 常 
数 ,y, 为 狄 拉克 箱 阵 ,y 是 夸克 化 学 势 年 阵 ,4 为 胶 子 
背景 场 ;第 二 项 是 介子 运动 学 项 ,Z,, 为 介子 波 函 数 
重 整 化 常数 ,是 介子 场 ;第 三 项 是 夸克 介子 相互 作 
用 项 ,hh 是 夸克 介子 之 间 的 Yukawa 相互 作用 强度 ， 
t= 1/2 (Is 地 :0),o 和 分 别 为 两 种 介子 的 场 :最 后 的 
两 项 分 别 为 介子 有 效 势 和 线性 破 缺 项 ,其 中 /p= 
1/2 内。 
在 上 面 所 介绍 的 低能 有 效 模型 的 计算 中 ,假设 
研究 的 体系 是 一 个 巨 正则 系 综 。 根 据 巨 正则 系 综 的 


p= Q[T, us] (59) 
如 果 是 要 计算 QCD 体系 的 状态 方程 ,可 以 计算 
压强 对 退 度 的 导数 ,从 而 得 到 体系 的 焙 等 热力 学 量 ， 
此 处 不 过 多 娄 述 ,下 面 主要 讨论 重子 数 涨 落 的 计算 。 
根据 巨 正则 系 综 的 性 质 ,粒子 数 的 涨 落 可 以 通 
过 求 压 强 对 化 学 势 的 各 阶 导数 得 到 。 在 不 考虑 电荷 
和 奇异 数 化 学 势 的 情况 下 , 重子 化 学 势 可 以 表示 为 
三 倍 的 夸克 化 学 势 心 = 34。 在 此 基础 上 ,可 以 得 到 
无 量 纲 的 重子 数 涨 落 ,也 就 是 广义 磁化 率 : 
ee o” 
an DOT a aD) 
WTR AE ARG Wy FEG BAS E E EA A 
Ou PAY oR SR a BE a EAA Td YP BY E 
子 数 涨 落 。 原 则 上 讲 , 只 要 拥有 足够 精确 的 压强 数 
据 , 就 可 以 对 任意 阶 的 重子 数 涨 落 进行 计算 ,但 是 受 
限于 计算 精度 ,计算 只 能 给 出 有 限 阶 的 重子 数 涨 落 ， 
目前 的 数值 方法 可 以 计算 到 第 十 阶 的 重子 数 涨 落 。 
计算 得 到 的 各 阶 重子 数 涨 落 还 不 能 直接 和 实验 
数据 作对 比 ,因为 实验 数据 直接 给 出 的 是 不 同 对 撞 
能 量 下 粒子 数 的 分 布 ,通过 其 数据 可 以 计算 出 该 分 
布下 的 各 阶 累积 量 , 而 累积 量 与 广义 磁化 率 之 间 还 
相差 一 个 体积 的 因子 ,如 下 所 示 : 


pa 


1 
xi = pr Na> 
1 
x2 =} <(6Nsy > 
VT 
i (61) 
X= yrs < (ONG) > 
wh =k (<(6N,) > 3 < (ôN) >?) 


YT 

这 里 给 出 了 从 第 一 到 第 四 阶 涨 落 与 广义 磁化 率 
的 关系 。 其 中 《 ) 代 表 统 计 平均 ,6N, =N,- 《Noy 
所 以 ,实际 上 不 可 以 直接 将 计算 得 到 的 广义 磁化 率 


基本 定义 ,我 们 的 体系 具有 确定 的 体积 和 温度 并 


与 实验 的 累积 量 进 行 对 比 。 为 了 消除 体积 因子 的 影 


可 以 与 外 热源 交换 粒子 和 能 量 。 这 也 就 是 为 什么 我 
们 的 标 度 依赖 的 有 效 作 用 量 中 引入 了 费 米 子 的 化 学 
势 。 在 巨 正则 系 综 的 框架 下 ,热力 学 量 的 计算 可 以 
归结 为 巨 热力 学 势 的 计算 。 对 于 处 于 热平衡 的 均匀 
系统 ,容易 得 到 巨 热力 学 势 密度 可 以 表示 成 红外 极 
Ek 一 oo 下 胶 子 势 . 介 子 有 效 势 和 线性 破 缺 项 的 和 : 

QUT, tg] = Veue(L,£) + Vn (p,L,L) -co (58) 

当 对 保 G 流 方程 求解 至 红外 后 就 可 以 得 到 式 
《58) 中 的 各 部 分 ,并 最 终 求 得 巨 热力 学 势 密度 。 此 
时 的 巨 热力 学 势 密度 依赖 于 温度 与 重子 化 学 势 ,也 
就 是 说 体系 的 热力 学 性 质 都 包含 在 其 中 。 在 此 基础 
上 根据 巨 正则 系 综 的 性 质 ,可 以 直接 得 到 所 研究 体 
系 的 压强 : 


响 , 可 以 将 广义 磁化 率 的 比值 与 累积 量 的 比值 进行 
比较 , 壁 如 : 偏 度 So = yy EE ko’ = y2 ly). 

62 ” 零 化 学 势 下 的 重子 数 涨 落 

图 7 给 出 了 零 化 学 势 下 重子 数 涨 落 的 比值 R2 = 
HEGRE, = rol? A0 RE, = Xe/2 随 光度 的 变化 行为 ， 
并 且 与 格 点 QCD 计算 的 结果 进行 比较 。 图 7 中 的 
红 带 给 出 了 了 促 G 低 能 有 效 理 论 计算 得 到 的 结果 , 绿 
带 和 蓝 带 给 出 的 是 格 点 HotQCD 合作 组 所 计算 的 结 
GRIT 1, 蓝 色 误差 棱 给 出 的 是 格 点 WB 合作 组 计 
算 的 结果 中 , 品 红色 虚线 给 出 的 是 强 子 共振 气体 模 
型 所 给 出 的 结果 "0。 这 里 保 G 结 果 中 的 误差 带 源 
自 于 对 温度 的 重 定 标 。 由 于 低能 有 效 理 论 与 第 一 性 
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图 7 FUZA PT BUI AY LE R? = x2/x3、RS = yë MRE = x2/x2 随 温度 的 变化 行为 ， 

并 且 与 格 点 QCD 计算 的 结果 进行 比较 外 

Fig.7 Ri, = 2/73, RE = xb /73 and Rs,= y#/%3 as functions of the temperature at vanishing baryon chemical 
potential (4, = 0) in comparison with lattice results” 


原理 QCD 的 计算 存在 温度 与 化 学 势 的 标 度 差异 ,所 加 ,所 以 ,有 限 化 学 势 的 计算 结果 对 低 碰撞 能 量 的 实 
以 引入 了 两 个 线性 重 定 标 系数 : 验 有 重要 的 意义 。 图 8 给 出 了 化 学 势 取 0~400 MeV 
TT = Cyl se 一 CU (62) ”情况 下 的 重子 数 涨 落 随 温度 的 依赖 关系 ,这 里 重子 

为 了 确定 这 两 个 系数 的 值 ,在 零 化 学 势 下 将 有 ”化 学 势 取 了 如 下 若干 数值 :ys 为 0 MeV、100 MeV, 

效 理论 的 结果 与 WB (Wuppertal-Budapest) 合 作 组 的 。 160 MeV、200 MeV、300 MeV 和 400 MeV。 涨 落 随 
格 点 QCD 结果 在 大 临 界 温度 处 进行 匹配 。 由 于 格 。” 着 化 学 势 升 高 的 行为 一 目 了 然 :各 阶 的 涨 落 都 在 较 
点 QCD 的 计算 存在 误差 ,这 就 导致 了 温度 重 定 标 系 ”高 的 化 学 势 下 产生 了 更 大 幅度 的 波动 。 此 外 , 随 着 
数 也 会 存在 相应 的 误差 。 对 于 化 学 势 重 定 标 系数 ， ”化 学 势 的 增 大 , 涨 落 在 高 于 相 变温 度 的 区 域 也 会 产 
这 可 以 通过 有 限 化 学 势 情况 下 匹配 相 边 界 的 曲率 来 。 生 一 个 新 的 峰 。 当 然 , 涨 落 的 阶 数 越 高 ,在 高 化 学 势 


确定 。 也 就 是 式 (63) 中 的 平方 项 的 曲率 系数 : 下 的 峰值 ( 谷 值 ) 会 变 得 越 大 ( 越 小 ) ,甚至 Re 在 化 学 
Tete) i e(a (63) B400 MeV 时 的 峰值 能 够 达到 2 000 左 右 。 可 以 
deya) (63) eee r 

ae. ep 注意 到 ,无 论 化 学 势 的 大 小 ,在 低温 区 域 总 是 
这 里 (po) 表 示 重 子 化 学 势 为 ,时 的 厦 临 界 po, -1, 这 是 由 于 在 低温 区 域 净重 子 数 的 分 布 趋 向 
温度 ,其 中 7. = Thua = 0)。 采 用 文献 [102] 中 的 相 。 于 Skellam 分 布 。 在 极端 高 温 区 域 ,由 于 QCD 渐 近 


边界 曲率 数值 x = 0.0153(18), 经 过 上 面 所 讨论 的 。 自由 的 性 质 ,夸克 和 胶 子 之 间 的 相互 作用 逐渐 趋 于 
匹配 过 程 , 可 以 得 到 两 个 重 定 标 系数 cr= 1.247(12) 。 消失 ,此 时 涨 落 的 取 值 也 逐渐 趋 于 Stefan-Boltzmann 
和 c,= 1.110(66), 这 两 个 系数 的 误差 是 保 G 有 限 温 ”极限 ,有 Ra(mn > 4) 一 0。 


度 有 限 化 学 势 下 计算 结果 的 误差 来 源 。 同时 也 容易 注意 到 , 当 化 学 势 增 大 时 ,计算 所 产 
从 图 7 可 以 看 到 ,经 过 重新 定 标 之 后 的 重子 数 。 生 的 误差 也 迅速 增加 ,可 以 看 到 ,在 高 阶 的 涨 落 中 高 


FN 


涨 落 从 四 阶 到 八 阶 都 与 格 点 QCD 的 计算 结果 符合 化 学 势 的 误差 几乎 填 满 了 曲线 所 在 的 区 域 。 这 说 明 


得 很 好 。 但 需要 注意 是 ,这 里 引入 线性 重 定 标 系 数 在 零 化 学 势 下 ,很 小 的 不 确定 性 会 导致 高 化 学 势 下 


的 方法 是 最 简单 ` 最 初步 地 将 低能 有 效 理论 与 第 一 ”明显 的 误差 ,这 导致 在 高 化 学 势 下 的 预言 具有 很 大 


性 原理 QCD 进行 联系 的 方法 。 之 所 以 采用 低能 有 ”的 不 确定 性 ,所 以 ,今后 的 一 个 重要 目标 是 尽量 减 小 
效 理论 来 计算 ,是 因为 其 简单 的 结构 能 够 相对 准确 a a a oer aga gett tp yp tate ty a 
ae 高 化 学 势 下 计算 的 误差 ,提高 理论 预言 的 可 靠 性 。 
地 刻画 整个 手 征 对 称 性 破 忽 与 恢复 的 过 程 。 当 然 ， 除了 温度 二 化 学 势 的 重 定 标 等 数 记 引 入 的 襄 关 
结构 简单 的 代价 是 缺失 了 高 能 标 处 规范 场 自由 度 所 a ae rs aR th do E RAIA Ea. 
带 来 的 动力 学 效应 ,所 以 更 合理 的 方法 是 在 绢 G。 。 从 前 面 给 出 的 专 克 介 子 模型 有 效 作用 量 可 以 看 到 ， 
wa R2: yee yee AX Yn sy -= ` 一 a J 

QCD 的 框架 下 ”计算 涨 落 等 观测 量 , 相 关 的 工作 1 目前 只 考虑 了 两 味 轻 专 克 以 及 标量 厦 标 介子 对 相 恋 
在 进行 中 。 在 下 节 中 可 以 很 自然 地 将 本 节 零 化 学 势 o OE REME ANETAR AEE 
的 计算 拓展 到 有 限 化 学 势 的 区 域 。 See veal ie 


T 


K Aaa 随 着 化 学 势 的 升 高 ,矢量 介子 以 及 di-quark 凝聚 等 
6.3 ”有 限 化 学 势 下 的 重子 数 涨 落 效应 都 会 对 体系 的 相 变 行为 产生 一 些 无 法 忽略 的 影 
在 重 离子 碰撞 实验 中 ,对 撞 能 量 的 降低 往往 标 。。 响 ,而 这 些 影响 显然 是 现 有 的 体系 所 没有 考虑 的 。 


志 着 体系 密度 的 增 大 ,也 就 对 应 着 体系 化 学 势 的 增 。 所 以 目前 在 高 化 学 势 区 域 进 行 的 计算 以 及 得 到 的 
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Fig.8 Re, RZ, and R$, as functions of the temperature with different values of ys in the range of 0~400 MeV”! 


些 结论 还 需要 谨慎 对 待 。 


所 以 目前 在 有 限 化 学 势 的 


计算 暂时 限制 在 实验 数据 所 覆盖 的 范围 内 ,也 就 是 
对 撞 能 量 在 Vsw = 7.7~200 GeV 的 范围 ,对 应 的 化 


学 势 约 在 心 = 20~400 MeV 的 范围 。 对 于 更 高 化 学 
势 区 域 涨 落 的 计算 ,需要 进一步 改进 低能 有 效 理论 
的 截断 ,在 其 中 引入 更 多 第 一 性 原理 QCD 计算 的 信 
息 , 相 关 的 工作 正在 进行 中 。 


相对 论 重 离子 碰撞 实验 中 核 核 对 撞 产 生 的 火 
球 ,经 过 热 化 形成 强 关 联 的 夸克 胶 子 等 离子 体 


(strongly intera 
sQGP 具 


近乎 完 美的 理想 流体 行为 rm 。 


cting Quark-Gluon Plasma, sQGP) ， 


SQGP 的 膨胀 及 其 演化 ,温度 不 断 降低 , 碎 列 成 各 种 


随 着 


各 样 的 强 子 , 从 而 实现 从 QGP 相 到 强 子 相 的 演化 ， 


但 是 ,此 时 强 子 


相互 作用 消失 ， 


之 间 还 存在 非 弹 性 相互 作用 , 强 子 的 


种 类 还 不 断 转 化 。 随 着 体系 的 进一步 演化 , 非 弹 性 


强 子 的 种 类 不 再 发 生变 化 ,这 一 时 刻 


所 对 应 的 过 程 有 一 个 专 有 名 称 ,被 称 为 化 学 冻结 。 
对 应 不 同 的 碰撞 能 量 , 化 学 冻结 时 刻 对 应 的 温度 和 


化 学 势 通常 是 不 同 的 。 显 然 ,实验 测量 的 涨 落 信号 


主要 受到 化 学 冻结 时 刻 体系 的 状态 所 影响 ,如 果 我 
们 进一步 忽略 非 平衡 演化 、 整 体重 子 数 守 恒 等 效 应 


的 潜在 影响 , 那 


么 ,只 要 我 们 在 不 同 对 撞 能 量 所 对 应 


的 化 学 冻结 的 温度 和 重子 化 学 势 上 计算 相应 的 涨 
落 , 就 可 以 和 实验 测量 进行 比较 。 这 里 化 学 冻结 的 
温度 和 重子 化 学 势 对 碰撞 能 量 的 依赖 关系 所 给 出 的 


1 线 就 是 化 学 冻结 曲线 。 


图 9 给 出 了 LHC 和 RHIC 等 实验 得 到 的 一 些 化 


学 冻结 温度 和 本 


E 子 化 学 势 的 实验 数据 点 以 及 若干 所 


使 用 的 化 学 冻结 曲线 。 蓝 色 五 角形 综合 
区 间 的 加 速 器 实验 所 得 到 的 化 学 冻结 实验 数据 ， 


各 个 能 量 


已 三 


红色 圆圈 代表 STAR 合作 组 得 到 的 化 学 冻结 数 


据 "” ,左下 角 的 小 图 放大 了 化 学 势 0-450 MeV 的 子 
区 间 。 图 9 中 的 蓝 色 虚 线 是 Andronic 等 拟 合 的 一 条 


化 学 冻结 曲线 " 


]。 它 的 参数 化 形式 可 以 写 为 ; 


a 
4 CF = 2 
"1 + 0.288 Vs 


(0) 
Ter 


T 


实验 数据 点 所 对 应 的 核子 对 质心 
中 的 红色 实 线 是 应 用 与 上 面相 同 


” 1+exp(2.60 - In(,/Smy )/0.45) 
其 中 :a = 1307.5 MeV, T& = 158.4 MeV; ,/ Syu 代表 


(64) 


ob E. 


能 量 ,GeV。 图 9 
的 参数 化 形式 到 


STAR 数据 点 上 进行 重新 拟 合 的 1 


1 线 。 这 里 用 了 


STAR 实验 所 给 出 的 全 部 7 个 数据 


点 来 重新 确定 式 


(64) 中 的 各 个 参数 。 对 于 STAR 合作 组 的 化 学 冻结 
数据 ,考虑 到 一 些 一 般 性 的 要 求 , 壁 如 化 学 冻结 温度 


随 化 学 冻结 化 学 势 的 增加 而 减 小 ， 
性 质 等 ,前 两 个 数据 点 以 及 最 后 一 


冻 出 曲线 的 凸凹 
个 数据 点 明显 有 


一 些 偏差 ,因此 ,更 合理 的 是 舍弃 这 些 不 合理 的 数据 


点 重新 进行 拟 合 , 这 样 就 得 到 了 图 9 中 第 三 条 化 学 


冻结 曲线 一 一 绿色 点 线 。 


图 10 给 出 了 温度 与 重子 化 学 


势 相 图 上 计算 得 


到 的 四 阶 与 二 阶 涨 落 的 比 R2, 此 外 ,图 10 还 给 出 了 


第 三 条 化 学 冻结 曲线 ,命名 为 STAR Fit I. WF 


200 


= 
= 120 


图 9 温度 与 化 学 势 平面 上 的 化 学 冻结 


线 以 及 Andronic 等 


和 STAR 合作 组 给 出 的 化 学 冻结 实验 数据 点 外 


Fig.9 Chemical freeze-out temperature a 


nd baryon chemical 


potential in the T - ys plane. The blue pentagons and red circles 
show the freeze-out data from Andronic et al and STAR 
experiment, respectively” 
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10 温度 与 化 学 势 平面 上 的 RS, 图 中 还 给 出 了 第 三 条 化 
学 冻结 曲线 ,命名 为 STAR Fit II 
Fig.10 Baryon number fluctuations R in the T - 4p plane as 
well as the freeze-out curve: STAR Fit II 


楚 地 看 到 , 随 着 化 学 势 的 增加 RS 逐渐 产生 更 加 剧烈 
的 波动 ,在 略微 高 于 相 变 温度 的 区 域 产生 了 负 值 的 
谷 , 这 同样 也 可 以 从 图 8 中 看 出 来 。 将 冻 出 曲线 画 
在 RS 的 相 图 上 ,可 以 看 到 , 冻 出 曲线 所 经 过 的 区 域 ， 
它 首先 靠近 了 红 蓝 交界 的 白色 区 域 然 后 又 弯 下 来 进 
入 了 深 红 色 区 域 。 这 样 的 行为 如 果 转 换 到 涨 落 随 化 
学 势 的 变化 图 中 就 能 够 看 到 一 个 随 化 学 势 的 非 单 调 
行为 。 

图 11 给 出 了 重子 数 涨 落 R2 、R2 和 R2 随 对 撞 能 
量 的 变化 行为 。 其 中 ,左边 两 条 误差 带 所 展示 的 结 


Hp /MeV 


SF PrP PG aE » 


Be RG-LEFT Freezeout: Andronic et al 
“=@™ {RG-LEFT Freezeout: STAR Fit I 
Å STAR Data (0%~5%) 


fRG-LEFT Freezeout: Andronic et al 
TAR Fit I 


{RG-LEFT Freezeout: S1 
STAR Data (0%~10%) 


N 9.9 Ga D dm N 
A AR YF YY ON 
V Syy / GeV 


果 对 应 着 我 们 在 图 9 中 所 采用 的 前 两 种 冻 出 曲线 的 
方案 。 灰 色 误 差 带 对 应 参数 化 的 冻 出 曲线 ,红色 误 
差 带 对 应 拟 合 STAR 数据 的 冻 出 曲线 ,误差 棒 代表 
实验 得 到 的 净 质 子 数 涨 落 的 数据 。 右 边 给 出 的 结 
是 使 用 了 第 三 条 化 学 冻结 曲线 :STAR Fit I。 在 网 
的 下 坐标 轴 给 出 了 实验 数据 所 对 应 的 对 撞 能 量 , 上 
坐标 轴 给 出 了 对 应 的 重子 化 学 势 。 对 于 RES 随 对 撞 
能 量变 化 的 行为 ,对 于 三 种 不 同 的 化 学 冻结 曲线 ,都 
发 现 随 着 碰撞 能 量 的 降低 ,R2 先 减少 后 增加 ,存在 
一 个 非 单调 性 ,这 与 实验 结果 是 一 致 的 。 这 种 非 单 
调 性 的 来 源 是 随 着 对 撞 能 量 的 降低 , 冻 出 化 学 势 越 
来 越 大 ,连续 过 渡 的 激烈 程度 越 来 越 强 , 这 样 就 导致 
了 RS 的 减少 ,明显 偏离 Skellam 分布 的 数值 ,但 是 随 
着 化 学 势 的 进一步 增加 ,化 学 冻结 曲线 开始 偏离 相 
边界 ,如 图 10 所 示 , 这 样 使 得 R2 进入 了 一 个 使 其 数 
明显 增加 的 区 域 。 对 于 六 阶 涨 落 ,理论 计算 发 现 
R2 在 很 大 的 能 量 区 间 是 负 值 , 这 与 实验 的 其 中 两 个 
数据 点 是 定性 符合 的 ,但 是 很 显然 需要 更 多 的 统计 
量 。 图 11 也 给 出 了 八 阶 涨 落 Re 的 理论 预言 ,但 是 
需要 注意 的 是 ,RS 的 误差 是 很 明显 的 ,特别 是 低能 
区 域 。 
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Fig.11 Baryon number fluctuations R$, RŽ, and R$, as functions of the collision energy obtained in QCD-assisted LEFT in 
comparison to the STAR data” 
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7 ”结语 


本 文 主 要 介绍 了 近年 来 保 G 在 QCD 相 变 及 其 


临界 现象 研究 中 的 一 些 进展 ,主要 侧重 于 与 QCD 的 


CEP 相关 的 


E 微 扰 临 界 指数 的 计算 、 重 子 数 的 涨 落 


等 。fRG 是 一 种 非 微 扰 连 续 场 论 的 理论 方法 , 它 提 
供 了 一 种 非常 适合 推导 重 整 化 群 流 方程 的 、 系 统 的 
HHE. FRG 既 可 以 用 来 做 系统 的 展开 ,譬如 梯 


度 展开 .e=4 
计算 。 


首先 ,介绍 了 应 用 保 G 研 究 dg 维 O(CV) 标 量 理 
论 临界 行为 的 一 些 重要 进展 。 利 用 保 G 很 容易 得 到 


- 4 展开 等 ,也 可 以 用 于 非 微 扰 场 论 的 


mn 


O(N ) 标 量 理论 非 微 扰 有 效 势 的 流 方程 。 这 个 非 微 


扰 有 效 势 的 流 方程 及 其 固定 点 方程 的 求解 ,对 于 顾 


究 二 级 相 变 的 临界 行为 至 关 重 要 。 一 方面 ,可 以 利 
用 有 效 势 泰 勒 展开 的 方法 对 其 展开 系数 的 流 方 程 逐 


一 求解 ,这 同 


Wilson 等 早年 发 展 的 a = 4 - dq 展开 的 


方法 在 本 质 上 是 一 致 的 。 在 名 G 的 框架 下 ,系统 地 


推导 且 重 复 了 临界 指数 se 展开 在 O(s) 阶 的 结果 , 原 


则 上 也 可 以 做 高 阶 的 计算 ,这 将 会 是 以 后 发 展 的 一 
个 方向 ,相应 的 研究 共有 很 重要 的 理论 意义 ,因为 
FRG 提供 了 一 个 可 以 简化 高 阶 计算 的 理论 框架 , 反 
过 来 ,se 展 开 的 微 扰 研究 也 可 以 促进 对 非 微 扰 固定 点 
方程 数学 结构 的 认识 。 值 得 注意 的 是 ,最 近 发 展 了 
一 种 不 依赖 于 微 扰 展开 ,直接 求解 非 微 扰 固定 点 方 


程 非 局 域 势 的 方法 ,该 方法 通过 解析 形式 考虑 了 有 


效 势 在 大 场 极限 情况 下 的 渐 近 行为 ,这 对 于 求解 空 


间 维 数 d 一 


2 的 临界 行为 至 关 重 要 ,也 是 保证 


MWH 定理 成 立 的 关键 。 这 种 非 局 域 势 直接 求解 的 


方法 具有 广阔 的 应 用 前 景 , 璧 如 动力 学 临界 指数 的 
计算 "5 、 多 固定 点 的 研究 "9 、Yang-Lee 边缘 奇 


(107-111 J Age 
点 等 。 


此 外 ,还 讨论 了 与 QCD 的 临界 终点 等 临界 现象 


密切 相关 的 


8 子 数 的 涨 落 , 这 是 相对 论 重 离子 碰撞 


实验 探寻 CE 


P 的 主要 观测 量 。 基 于 保 G 的 理论 计 


算 ,探讨 了 实验 观测 到 的 净 质 子 数 分 布 的 峰 度 ,也 就 


是 四 阶 矩 随 碰 撞 能 量 非 单调 依赖 关系 的 可 能 原因 : 
一 方面 , 随 着 碰撞 能 量 的 减 小 ,化 学 冻结 重子 化 学 势 


不 断 增加 ,导致 QCD 相 图 上 的 手 征 连 续 过 渡 越 来 越 


剧烈 , 涨 落 的 


幅度 越 来 越 大 ; 男 一 方面 ,化 学 冻结 曲 


线 在 高 化 学 势 的 区 域 越 来 越 往 低温 的 方向 偏离 相 边 
界 曲 线 。 为 了 最 终 确定 CEP 在 QCD 相 图 上 的 位 置 ， 


4 


不 仅 需 要 不 断 降低 高 化 学 势 区 域 理 论 计算 的 误差 ， 
比如 CEP 位 置 的 理论 预言 、 涨 落 的 计算 等 等 ,而 且 


也 需要 实验 提供 更 多 的 ,特别 是 高 化 学 势 区 域 的 化 


学 冻结 数据 。 


致谢 “感谢 北京 大 学 刘 玉 侈 教授 、 复 旦 大 学 马 余 刚 
院士 .华中 师范 大 学 罗 晓 峰 教 授 的 邀请 和 鼓励 ,使 得 
这 篇 综述 得 以 顺利 完成 ! 此 外 ,感谢 我 们 多 年 来 的 
合作 者 :Jens Braun, Jan M. Pawlowski, Fabian 
Rennecke, Bernd-Jochen Schaefer, Nicolas Wink, 陈 
AA Ba WEA, FRE, SHR DA. 

作者 贡献 声明 妖 诗 计算 并 撰写 文章 的 重子 数 涨 落 
部 分 ,修改 文章 ; 谈 阳 阳 计 算 并 撰写 文章 临界 指数 部 
分 ,文章 排版 ; 付 伟 杰 撰写 文章 引言 以 及 理论 介绍 部 
分 ,指导 修改 文章 。 


1 Luo X F, Wang Q, Xu N, et al. Properties of QCD matter 
at high baryon density[M]. Singapore: Springer Nature 
Singapore, 2022. DOI: 10.1007/978-98 1-19-4441-3. 

RB BE, REE. QCD 临界 点 附近 的 动力 学 临界 涨 落 
[J]. 核 技 术 , 2023, 46(4): 040004. DOI: 10.11889/j.0253- 
3219.2023.hjs.46.040004. 

WU Shanjin, SONG Huichao. Critical dynamical 


N 


fluctuations near the QCD critical point[J]. Nuclear 
Techniques, 2023, 46(4): 040004. DOI: 10.11889/j.0253- 
3219.2023.hjs.46.040004. 

3 UES, 黄梅. QCD Imi HRA A SH BHT]. 核 技术 ， 
2023, 46(4): 040005. DOI: 10.11889/j. 0253-3219.2023. 
hjs.46.040005. 

XU Kun, HUANG Mei. QCD critical end point and 


baryon number fluctuation[J]. Nuclear Techniques, 2023, 
46(4): 040005. DOI: 10.11889/. 0253-3219.2023. hjs. 
46.040005. 

4 SR, 张 定 伟 , 罗 晓 峰 . 相对 论 重 离子 碰撞 中 QCD 相 图 
的 实验 研究 四. 核 技 术 , 2023, 46(4): 040001. DOI: 
10.11889/j.0253-3219.2023.hjs.46.040001. 

ZHANG Yu, ZHANG Dingwei, LUO Xiaofeng. 


Experimental study of the QCD phase diagram in 
relativistic heavy-ion collisions[J]. Nuclear Techniques, 
2023, 46(4): 040001. DOI: 10.11889/j. 0253-3219.2023. 
hjs.46.040001. 

5 ， 杜 轶 伦 , 李 程 明 , 史 潮 , 等 . 基于 有 效 场 论 的 QCD HAE 
Ft [J]. 核 技术 ，2023，46(4): 040009. DOT: 10.11889/ 
j.0253-3219.2023.hjs.46.040009. 

DU Yilun, LI Chengming, SHI Chao, et al. Review of 


QCD phase diagram analysis using effective field theories 
[J]. Nuclear Techniques, 2023, 46(4): 040009. DOI: 
10.11889/j.0253-3219.2023.hjs.46.040009. 

6 Aoki Y, Endrodi G, Fodor Z, et al. The order of the 


040002-13 


202306.00144v1 


chinaXiv 


核 技 术 


ChinaXiv 合 作 期 刊 


2023, 46: 040002 


10 


11 


12 


13 


14 


15 


16 


17 


quantum chromodynamics transition predicted by the 
standard model of particle physics[J]. Nature, 2006, 443 
(7112): 675 - 678. DOI: 10.1038/nature05120. 

Andronic A, Braun-Munzinger P, Redlich K, et al. 
Decoding the phase structure of QCD via particle 
production at high energy[J]. Nature, 2018, 561(7723): 
321 - 330. DOI: 10.1038/s41586-018-0491-6. 

Borsanyi S, fodor Z, Giordano M, et al. Equation of state 
of a hot-and-dense quark gluon plasma: lattice 
simulations at real u, vs. extrapolations[J]. 2022. DOI: 
10.48550/arXiv.2208.05398. 

Borsanyi S, Fodor Z, Guenther J N, et al. Lattice QCD 
equation of state at finite chemical potential from an 
alternative expansion scheme[J]. Physical Review Letters, 
2021, 126(23): 232001. DOI: 10.1103/physrevlett. 126. 
232001. 

Bollweg D, Goswami J, Kaczmarek O, et al. Taylor 
expansions and Padé approximants for cumulants of 
conserved charge fluctuations at nonvanishing chemical 
potentials[J]. Physical Review D, 2022, 105(7): 074511. 
DOI: 10.1103/physrevd.105.074511. 

Fischer C S. QCD at finite temperature and chemical 
potential from Dyson-Schwinger equations[J]. Progress in 
Particle and Nuclear Physics, 2019, 105: 1 - 60. DOI: 10. 
1016/j.ppnp.2019.01.002. 

Fu W J, Pawlowski J M, Rennecke F. QCD phase 
structure at finite temperature and density[J]. Physical 
Review D, 2020, 101(5): 054032. DOI: 10.1103/ 
physrevd.101.054032. 

Gao F, Pawlowski J M. Chiral phase structure and critical 
end point in QCDIJ]. Physics Letters B, 2021, 820: 
136584. DOI: 10.1016/j.physletb.202 1.136584. 

Gunkel P J, Fischer C S. Locating the critical endpoint of 
QCD: Mesonic backcoupling effects[J]. Physical Review 
D, 2021, 104(5): 054022. DOI: 10.1103/physrevd. 104. 
054022. 

Fu W J. QCD at finite temperature and density within the 
fRG approach: Communications in 
Theoretical Physics, 2022, 74(9): 097304. DOI: 10.1088/ 
1572-9494/ac86be. 

Schaefer B J, Wambach J. The phase diagram of the 
quark-meson model[J]. Nuclear Physics A, 2005, 757(3 - 
4): 479 - 492. DOI: 10.1016/j.nuclphysa.2005.04.012. 


an overview|[J]. 


Fukushima K, Skokov V. Polyakov loop modeling for hot 
QCDIJ]. Progress in Particle and Nuclear Physics, 2017, 
96: 154 - 199. DOI: 10.1016/j.ppnp.2017.05.002. 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


040002-14 


Chen Y R, Wen R, Fu W J. Critical behaviors of the O(4) 
and Z(2) symmetries in the QCD phase diagram[J]. 
Physical Review D, 2021, 104(5): 054009. DOI: 10.1103/ 
physrevd.104.054009. 

Grossi E, Ihssen F J, Pawlowski J M, et al. Shocks and 
quark-meson scatterings at large density[J]. Physical 
Review D, 2021, 104: 016028. DOI: 10.1103/physrevd. 
104.016028. 

Otto K, Busch C, Schaefer B J. Regulator scheme 
dependence of the chiral phase transition at high densities 
[J]. Physical Review D, 2022, 106(9): 094018. DOI: 
10.1103/physrevd.106.094018. 

Ma S G. Modem theory of critical phenomena[M]. 
Cambridge, Mass: Perseus Publishing, 2000. 

Wilson K G. Renormalization group and critical 
phenomena. I. renormalization group and the kadanoff 
scaling picture[J]. Physical Review B, 1971, 4(9): 3174 - 
3183. DOI: 10.1103/physrevb.4.3174. 

Wilson K G. Renormalization group and critical 
phenomena. II. phase-space cell analysis of critical 
behavior[J]. Physical Review B, 1971, 4(9): 3184 - 3205. 
DOI: 10.1103/physrevb.4.3184. 

Wilson K G, Fisher M E. Critical exponents in 3.99 
dimensions[J]. Physical Review Letters, 1972, 28(4): 
240 - 243. DOI: 10.1103/physrevlett.28.240. 

Wilson K. The renormalization group and the gs expansion 
[J]. Physics Reports, 1974, 12(2): 75 - 199. DOI: 10.1016/ 
0370-1573(74)90023-4. 

Luo X F, Xu N. Search for the QCD critical point with 
fluctuations of conserved quantities in relativistic heavy- 
ion collisions at RHIC: an overview[J]. Nuclear Science 
and Techniques, 2017, 28(8): 112. DOI: 10.1007/s41365- 
017-0257-0. 

Bzdak A, Esumi S, Koch V, et al. Mapping the phases of 
quantum chromodynamics with beam energy Scan[J]. 
Physics Reports, 2020, 853: 1 - 87. DOI: 10.1016/. 
physrep.2020.01.005. 

Friman B, Höhne C, Knoll J, et al. The CBM physics 
book: 
experiments[M]. 
Heidelberg, 2011. DOI: 10.1007/978-3-642-13293-3. 
Agakishiev G, Balanda A, Bannier B, et al. The high- 


compressed baryonic matter in laboratory 


Berlin, Heidelberg: Springer Berlin 


acceptance dielectron spectrometer HADES[J]. European 
Physical Journal, 2009, A41: 243 - 277. DOI: 10.1140/ 
epja/i2009-10807-5. 

Abgrall N, Andreeva O, Aduszkiewicz A, et al. NA61/ 


202306.00144v1 


chinaXiv 


临界 现象 与 泛 函 重 整 化 群 


ChinaXiv 合 作 期 刊 


31 


32 


33 


34 


35 


36 


37 


38 


39 


40 


41 


42 


SHINE facility at the CERN SPS: beams and detector 
system[J]. JINST, 2014, 9: P06005. DOI: 10.1088/1748- 
0221/9/06/P06005. 

Sorin A, Kekelidze V, Kovalenko A, et al. Heavy-ion 
program at NICA/MPD at JINR[J]. Nuclear Physics A, 
2011, 855(1): 510 - 513. DOI: 10.1016/.nuclphysa.2011. 
02.118. 

Blaschke D, Aichelin J, Bratkovskaya E, et al. Topical 
issue on exploring strongly interacting matter at high 
densities - NICA white paper[J]. The European Physical 
Journal A, 2016, 52(8): 267. DOL 10.1140/epja/i2016- 
16267-x. 

Dainese A, Diehl M, Di Nezza P, et al. Physics beyond 
colliders: QCD working group report[EB/OL]. 2019: 
arXiv: 1901.04482. https://arxiv.org/abs/1901.04482. 

Yang J C, Xia J W, Xiao G Q, et al. High Intensity heavy 
ion Accelerator Facility (HIAF) in China[J/OL]. Nuclear 
Instruments & Methods, 2013, B317: 263 - 265. DOI: 
10.1016/j.nimb. 2013.08.046. 

Lü L M, Yi H, Xiao Z G, et al. Conceptual design of the 
HIRFL-CSR external-target experiment[J]. Science China 
Physics, Mechanics & Astronomy, 2017, 60(1): 012021. 
DOI: 10.1007/s11433-016-0342-x. 

Sakaguchi T. Study of high baryon density QCD matter at 
J-PARC-HI[J]. Nuclear Physics A, 2017, 967: 896 - 899. 
DOI: 10.1016/j.nuclphysa.2017.05.081. 

Sako H. Studies of extremely dense matter in heavy-ion 
collisions at J-PARC[J]. Nuclear Physics A, 2019, 982: 
959 - 962. DOL: 10.1016/j.nuclphysa.2018.11.027. 
Stephanov M A, Rajagopal K, Shuryak E V. Event-by- 
event fluctuations in heavy ion collisions and the QCD 
critical point[J]. Physical Review D, 1999, 60(11): 
114028. DOI: 10.1103/physrevd.60.114028. 

Stephanov M A. Non-Gaussian fluctuations near the QCD 
critical point[J]. Physical Review Letters, 2009, 102(3): 
032301. DOI: 10.1103/physrevlett.102.032301. 

Stephanov M A. Sign of kurtosis near the QCD critical 
point[J]. Physical Review Letters, 2011, 107(5): 052301. 
DOI: 10.1103/physrevlett.107.052301. 

Adamczyk L, Adkins J K, Agakishiev G, et al. Energy 
dependence of moments of net-proton multiplicity 
distributions at RHIC[J]. Physical Review Letters, 2014, 
112: 032302. DOI: 10.1103/PhysRevLett.112.032302. 
Adamczyk L, Adkins J K, Agakishiev G, et al. Beam 
energy dependence of moments of the net-charge 


multiplicity distributions in Au+Au collisions at RHIC[J]. 


43 


44 


45 


46 


47 


48 


49 


50 


51 


52 


53 


040002-15 


Physical Review Letters, 2014, 113: 092301. DOI: 
10.1103/PhysRevLett.113.092301. 

Luo X. Energy dependence of moments of net-proton and 
net-charge multiplicity distributions at STAR[J]. PoS, 
2015, CPOD2014: 019. 

Adamezyk L, OTHERS. Collision energy dependence of 
moments of net-kaon multiplicity distributions at RHIC 
[J]. Physics Letters, 2018, B785: 551 - 560. DOI: 10. 
1016/j.physletb.2018.07.066. 

Adam J, Adamczyk L, Adams J R, et al. Collision-energy 
dependence of second-order off-diagonal and diagonal 
cumulants of net-charge, net-proton, and net-kaon 
multiplicity distributions in Au+Au collisions[J]. Physical 
Review C, 2019, 100(1): 014902. DOI: 10.1103/ 
PhysRevC.100.014902. 

Adam J, Adamczyk L, Adams J R, et al. Nonmonotonic 
energy dependence of net-proton number fluctuations[J]. 
Physical Review Letters, 2021, 126(9): 092301. DOI: 
10.1103/physrevlett.126.092301. 

Abdallah M, Adam J, Adamczyk L, et al. Measurement of 
the sixth-order cumulant of net-proton multiplicity 
distributions in Au+Au collisions at J Syn= 27, 54.4, and 
200 GeV at RHIC[J]. Physical Review Letters, 2021, 127 
(26): 262301. DOI: 10.1103/PhysRevLett.127.262301. 
Abdallah M S, Aboona B E, Adam J, et al. Measurements 
of proton high-order cumulants in ,/S,,=3 GeV AutAu 
collisions and implications for the QCD critical point[J]. 
Physical Review Letters, 2022, 128(20): 202303. DOI: 
10.1103/physrevlett.128.202303. 

Karsch F. Critical behavior and net-charge fluctuations 
from lattice QCD[J]. PoS, 2019, CORFU2018: 163. DOI: 
10.22323/1.347.0163. 

Wetterich C. Exact evolution equation for the effective 
potential[J]. Physics Letters B, 1993, 301(1): 90 - 94. 
DOI: 10.1016/0370-2693(93)90726-x. 

Berges J, Tetradis N, Wetterich C. Non-perturbative 
renormalization flow in quantum field theory and 
statistical physics[J]. Physics Reports, 2002, 363(4 - 6): 
223 - 386. DOI: 10.1016/s0370-1573(01)00098-9. 
Pawlowski J M. Aspects of the functional renormalisation 
group[J]. Annals of Physics, 2007, 322(12): 2831 - 2915. 
DOI: 10.1016/j.aop.2007.01.007. 

Schaefer B J, Wambach J. Renormalization group 
approach towards the QCD phase diagram[J]. Physics of 
Particles and Nuclei, 2008, 39(7): 1025 - 1032. DOI: 


202306.00144v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


核 技 术 2023,46: 040002 
10.1134/S1063779608070083. DOI: 10.1016/s0370-2693(00)00748-6. 
54 Gies H. Introduction to the functional RG and 67 Litim D F. Optimized renormalization group flows[J]. 


J3 


56 


57 


58 


59 


60 


61 


62 


63 


64 


65 


66 


applications to gauge theories[J]. Lecture Notes in 
Physics, 2012, 852: 287 - 348. DOI: 10.1007/978-3-642- 
27320-9_6. 

Rosten O J. Fundamentals of the exact renormalization 
group[J]. Physics Reports, 2012, 511(4): 177 - 272. DOI: 
10.1016/j.physrep.2011.12.003. 

Braun J. Fermion interactions and universal behavior in 
strongly interacting theories[J]. Journal of Physics G: 
Nuclear and Particle Physics, 2012, 39(3): 033001. DOI: 
10.1088/0954-3899/39/3/033001. 

Pawlowski J M. Equation of state and phase diagram of 
strongly interacting matter[J]. Nuclear Physics A, 2014, 
931: 113 - 124. DOI: 10.1016/j.nuclphysa.2014.09.074. 
Dupuis N, Canet L, Eichhom A, et al. The 
nonperturbative functional renormalization group and its 
applications[J]. Physics Reports, 2021, 910: 1 - 114. DOI: 
10.1016/j.physrep.2021.01.001. 

Braun J, Fister L, Pawlowski J M, et al. From quarks and 
gluons Chiral 
dynamical QCD[J]. Physical Review D, 2016, 94(3): 
034016. DOI: 10.1103/physrevd.94.034016. 

Mitter M, Pawlowski J M, Strodthoff N. Chiral symmetry 


to hadrons: symmetry breaking in 


breaking in continuum QCD{J]. Physical Review D, 
2015, 91(5): 054035. DOI: 10.1103/physrevd.91.054035. 
Rennecke F. Vacuum structure of vector mesons in QCD 
[J]. Physical Review D, 2015, 92(7): 076012. DOI: 10. 
1103/physrevd.92.076012. 

Cyrol A K, Fister L, Mitter M, et al. Landau gauge Yang- 
Mills correlation functions[J]. Physical Review D, 2016, 
94(5): 054005. DOI: 10.1103/physrevd.94.054005. 

Cyrol A K, Mitter M, Pawlowski J M, ef al. 
Nonperturbative finite-temperature Yang-Mills theory[J]. 
Physical Review D, 2018, 97(5): 054015. DOI: 10.1103/ 
physrevd.97.054015. 

Cyrol A K, Mitter M, Pawlowski J M, et al. 
Nonperturbative quark, gluon, and meson correlators of 
unquenched QCD[J]. Physical Review D, 2018, 97(5): 
054006. DOI: 10.1103/physrevd.97.054006. 

Braun J, Fu W J, Pawlowski J M, et al. Chiral 
susceptibility in (2+1) -flavor QCD[J]. Physical Review 
D, 2020, 102(5): 056010. DOI: 10.1103/physrevd. 102. 
056010. 

Litim D F. Optimisation of the exact renormalisation 


group[J]. Physics Letters B, 2000, 486(1 - 2): 92 - 99. 


68 


69 


70 


71 


72 


73 


74 


75 


76 


77 


78 


040002-16 


Physical Review D, 2001, 64(10): 105007. DOI: 10.1103/ 
physrevd.64.105007. 
Litim D F. Critical exponents from optimised 
renormalisation group flows[J]. Nuclear Physics B, 2002, 
631(1 - 2): 128 - 158. DOI: 10.1016/s0550-3213(02) 
00186-4. 

Mermin N D, Wagner H. Absence of ferromagnetism or 
antiferromagnetism in one- or two-dimensional isotropic 
Heisenberg models[J]. Physical Review Letters, 1966, 17 
(22): 1133 - 1136. DOI: 10.1103/physrevlett.17.1133. 
Hohenberg P C. Existence of long-range order in one and 
two dimensions[J]. Physical Review, 1967, 158(2): 383 - 
386. DOI: 10.1103/physrev. 158.383. 

Coleman S R. There are no Goldstone bosons in two 
dimensions[J]. Communications in Mathematical Physics, 
1973, 31(4): 259 - 264. DOI: 10.1007/BF01646487. 

Tan Y Y, Huang C, Chen Y R, et al. Criticality of the O 
(N) universality via global solutions to nonperturbative 
fixed-point equations[EB/OL]. 2022: arXiv: 2211.10249. 
https://arxiv.org/abs/2211.10249. 

Bazavov A, Ding H T, Hegde P, et al. Freeze-out 
from QCD 
thermodynamics[J]. Physical Review Letters, 2012, 109 


(19): 192302. DOI: 10.1103/physrevlett.109.192302. 


conditions in heavy ion collisions 


Borsanyi S, Fodor Z, Katz S D, et al. Freeze-out 
parameters: lattice meets experiment[J]. Physical Review 
Letters, 2013, 111(6): 062005. DOI: 10.1103/physrevlett. 
111.062005. 

Borsanyi S, Fodor Z, Katz S D, et al. Freeze-out 
parameters from electric charge and baryon number 
fluctuations: is there consistency? [J]. Physical Review 
Letters, 2014, 113(5): 052301. DOI:  10.1103/ 
PhysRevLett.113.052301. 

Bazavov A, Ding H T, Hegde P, et al. Skewness and 
kurtosis of net baryon-number distributions at small 
values of the baryon chemical potential[J]. Physical 
Review D, 2017, 96(7): 074510. DOI: 10.1103/physrevd. 
96.0745 10. 

Bazavov A, Ding H T, Hegde P, et al. QCD equation of 
state to O(u,”) from lattice QCD[J]. Physical Review D, 
2017, 95(5): 054504. DOI: 10.1103/physrevd.95.054504. 
Bazavov A, Bollweg D, Ding H T, et al. Skewness, 
kurtosis, and the fifth and sixth order cumulants of net 


baryon-number distributions from lattice QCD confront 


202306.00144v1 


chinaXiv 


临界 现象 与 泛 函 重 整 化 群 


ChinaXiv 合 作 期 刊 


79 


80 


81 


82 


83 


84 


85 


86 


87 


88 


89 


high-statistics STAR data[J]. Physical Review D, 2020, 
101(7): 074502. DOI: 10.1103/physrevd.101.074502. 
Borsanyi S, Fodor Z, Guenther J N, et al. Higher order 
fluctuations and correlations of conserved charges from 
lattice QCD[J]. Journal of High Energy Physics, 2018, 
2018(10): 205. DOI: 10.1007/JHEP10(2018)205. 

Skokov V, Stokić B, Friman B, et al. Meson fluctuations 
and thermodynamics of the Polyakov-loop-extended 
quark-meson model[J]. Physical Review C, 2010, 82: 
015206. DOI: 10.1103/physrevc.82.015206. 

Skokov V, Friman B, Redlich K. Quark number 
fluctuations in the Polyakov loop-extended quark-meson 
model at finite baryon density[J]. Physical Review C, 
2011, 83(5): 054904. DOI: 10.1103/physrevc.83.054904. 
Friman B, Karsch F, Redlich K, et al. Fluctuations as 
probe of the QCD phase transition and freeze-out in 
heavy ion collisions at LHC and RHIC[J]. The European 
Physical Journal C, 2011, 71(7): 1694. DOI: 10.1140/epjc/ 
s10052-011-1694-2. 

Morita K, Friman B, Redlich K. Criticality of the net- 
baryon number probability distribution at finite density[J]. 
Physics Letters B, 2015, 741: 178 - 183. DOI: 10.1016/j. 
physletb.2014. 12.037. 

Fu W J, Pawlowski J M. Relevance of matter and glue 
dynamics for baryon number fluctuations[J]. Physical 
Review D, 2015, 92(11): 116006. DOI: 10.1103/physrevd. 
92.116006. 

Fu W J, Pawlowski J M. Correlating the skewness and 
kurtosis of baryon number distributions[J]. Physical 
Review D, 2016, 93(9): 091501. DOI: 10.1103/physrevd. 
93.091501. 

Fu W J, Pawlowski J M, Rennecke F, et al. Baryon 
number fluctuations at finite temperature and density[J]. 
Physical Review D, 2016, 94(11): 116020. DOI: 10.1103/ 
physrevd.94.116020. 

Almasi G A, Friman B, Redlich K. Baryon number 
fluctuations in chiral effective models and their 
phenomenological implications[J]. Physical Review D, 
2017, 96(1): 014027. DOI: 10.1103/physrevd.96.014027. 
Fu W J, Pawlowski J M, Rennecke F. Strangeness 
neutrality and QCD thermodynamics[J]. SciPost Physics 
Core, 2020, 2: 2. DOI: 10.21468/scipostphyscore.2. 1.002. 
Fu W J, Pawlowski J M, Rennecke F. Strangeness 
neutrality and 
Physical Review D, 2019, 100(11): 111501. DOI: 10.1103/ 
physrevd.100.111501. 


baryon-strangeness _correlations[J]. 


90 


91 


92 


93 


94 


95 


96 


97 


98 


99 


100 


101 


040002-17 


Wen R, Huang C, Fu W J. Baryon number fluctuations in 
the 2+1 flavor low energy effective model[J]. Physical 
Review D, 2019, 99(9): 094019. DOI: 10.1103/physrevd. 
99.094019. 

Wen R, Fu W J. Correlations of conserved charges and 
QCD phase structure[J]. Chinese Physics C, 2021, 45(4): 
044112. DOI: 10.1088/1674-1137/abe199. 

Fu W J, Luo X F, Pawlowski J M, et al. Hyper-order 
baryon number fluctuations at finite temperature and 
density[J]. Physical Review D, 2021, 104(9): 094047. 
DOI: 10.1103/physrevd. 104.094047. 

Fu W J, Liu Y X, Wu Y L. Fluctuations and correlations 
of conserved charges in QCD at finite temperature with 
effective models[J]. Physical Review D, 2010, 81: 
014028. DOI: 10.1103/physrevd.81.014028. 

Fu W J, Wu Y L. Fluctuations and correlations of 
conserved charges near the QCD critical point[J]. 
Physical Review D, 2010, 82(7): 074013. DOI: 10.1103/ 
physrevd.82.074013. 

Karsch F, Schaefer B J, Wagner M, et al. Towards finite 
density QCD with Taylor expansions[J]. Physics Letters 
B, 2011, 698(3): 256 - 264. DOI: 10.1016/j. physletb. 
2011.03.013. 

Schaefer B J, Wagner M. QCD critical region and higher 
moments for three-flavor models[J]. Physical Review D, 
2012, 85(3): 034027. DOI: 10.1103/physrevd.85.034027. 
Li Z B, Xu K, Wang X Y, et al. The kurtosis of net baryon 
number fluctuations from a realistic Polyakov-Nambu- 
Jona-Lasinio model along the experimental freeze-out line 
[J]. The European Physical Journal C, 2019, 79(3): 245. 
DOI: 10.1140/epjc/s10052-019-6703-x. 

Xin X Y, Qin S X, Liu Y X. Quark number fluctuations at 
finite temperature and finite chemical potential via the 
Dyson-Schwinger equation approach[J]. Physical Review 
D, 2014, 90(7): 076006. DOI: 10.1103/physrevd. 90. 
076006. 

Isserstedt P, Buballa M, Fischer C S, et al. Baryon 
number fluctuations in the QCD phase diagram from 
Dyson-Schwinger equations[J]. Physical Review D, 2019, 
100(7): 074011. DOI: 10.1103/physrevd.100.074011. 
Bazavov A, Ding H T, Hegde P, et al. Skewness and 
kurtosis of net baryon-number distributions at small 
values of the baryon chemical potential[J]. Physical 
Review D, 2017, 96(7): 074510. DOI: 10.1103/physrevd. 
96.0745 10. 

Braun-Munzinger P, Redlich K, Stachel J. Particle 


202306.00144v1 


chinaXiv 


核 技 术 


ChinaXiv 合 作 期 刊 


2023, 46: 040002 


102 


103 


104 


105 


106 


production in heavy ion collisions[M]. Quark-Gluon 
Plasma 3. Singapore: World Scientific, 2004: 491 - 599. 
DOI: 10.1142/9789812795533_0008. 

Borsanyi S, Fodor Z, Guenther J N, et al. QCD crossover 
at finite chemical potential from lattice simulations[J]. 
Physical Review Letters, 2020, 125(5): 052001. DOI: 10. 
1103/PhysRevLett.125.052001. 

Hwa R C, Wang X N. Quark-gluon plasma 3[M]. 
Singapore: World Scientific, 2004. DOI: 10.1142/5029. 
Adamczyk L, Adkins J K, Agakishiev G, et al. Bulk 
properties of the medium produced in relativistic heavy- 
ion collisions from the beam energy scan program{[J]. 
Physical Review C, 2017, 96(4): 044904. DOI: 10.1103/ 
PhysRevC.96.044904. 

Tan Y Y, Chen Y R, Fu W J. Real-time dynamics of the O 
(4) scalar theory within the fRG approach[J]. SciPost 
Physics, 2022, 12: 26. DOI: 10.21468/scipostphys. 12. 
1.026. 

Yabunaka S, Delamotte B. Surprises in O(N) models: 
nonperturbative fixed points, large N limits, and 


multicriticality[J]. Physical Review Letters, 2017, 119 


107 


108 


109 


110 


111 


040002-18 


(19): 191602. DOI: 10.1103/PhysRevLett.119.191602. 
Stephanov M A. QCD critical point and complex 
chemical potential singularities[J]. Physical Review D, 
2006, 73(9): 094508. DOI: 10.1103/physrevd.73.094508. 
Mukherjee S, Skokov V. Universality driven analytic 
structure of the QCD crossover: radius of convergence in 
the baryon chemical potential[J]. Physical Review D, 
2021, 103(7): L071501. DOI: 10.1103/physrevd. 103. 
1071501. 

Connelly A, Johnson G, Rennecke F, et al. Universal 
location of the Yang-lee edge singularity in O(N) theories 
[J]. Physical Review Letters, 2020, 125(19): 191602. 
DOI: 10.1103/PhysRevLett. 125.191602. 

Rennecke F, Skokov V. Universal location of Yang-lee 
edge singularity for a one-component field theory in 1<D< 
4[J]. Annals of Physics, 2022, 444: 169010. DOI: 
10.1016/j.aop.2022.169010. 

Thssen F, Pawlowski J M. Functional flows for complex 
effective actions[EB/OL]. 2022: arXiv: 2207.10057. https: 
//arxiv.org/abs/2207. 10057. 


